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On the bispectrum of COBE and WMAP 
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ABSTRACT 

The COBE-DMR 4-year maps displayed a strong non-Gaussian signal in the "inter- 
scale" components of the bispectrum: their observed values did not display the scatter 
expected from Gaussian maps. We re-examine this and other suggested non-Gaussian 
features in the light of WMAP. We find that they all disappear. Given that it was 
proved that COBE-DMR high noise levels and documented systematics could at most 
dilute the observed non-Gaussian features, we conclude that this dataset must have 
contained non-negligible undocumented systematic errors. It turns out that the culprit 
is a combination of QuadCube pixelization and data collected during the "eclipse 
season" . 

Key words: cosmic microwave background - Gaussianity tests. 



1 INTRODUCTION 

The possibility of non-Gaussianity in the COBE-DMR 
4 year maps led to a rat her protracted story. Using 
"singl e-l" bispectrum analysis. iFerreira. Maeueiio fc Gorskil 
( 1998) found strong evidence for non-Gaussianity in 
the anisotropy of the Cosmic Microwave Background 
(CMB) temperature. T his detection was followed by 
simila r clai ms bv | N ovikov 1 Feldman and S handarir] 
Jl99Sft and iPando. Valls-Gabaud fc Fand Jl99^ and 
caused con siderabl e con s ternat ion among theorists (see 
iKamionkowski and Jaffel <ll998l) for a discussion). Fur- 
ther work, however, showed that claims other than 
those based on the bispectrum coul d not be reproduced 
jMukheriee. Hobson. Las enby| | 2000ft. The bispe ctrum 
claims wer e confirmed bylBromlex^^Sksiark| dl999T) . 

Later iBandav. Zaroubi^^G^rskiri^99l) cast serious 
doubts upon the cosmological origin of the observed sig- 
nal. A systematic was identified which removed the observed 
"single-^"' bispectrum signal - the so-called eclipse season 
data, which should never have been used. However, when an 
extension to "inter-^" bispectrum components was sought, a 
new non -Gaussian signa l was found for inter-^ separations of 
Al = 1 iMagueiidl20 0C0 . Specifically, their observed values 
were found to concentrate uncannily close to zero instead of 
displaying the scatter expected from Gaussian maps. 
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(Maguciio 2000). It was also found that the observed sig- 
nal did not extend to hig her inter-^ separations Al > 1 
JSandvikfc MagueiidE)0lT> . 

What shall we ma ke of these claims in the light of 
the recent observations feennett fc a j| |2003fr by the Wilkin- 
son Microwave Anisotropy Probe (WMAP)? In this paper 
we show that they don't survive the new data, which dis- 
plays consistency with Gaussianity on the angular scales 
probed by COBE. This might seem at first suprising, given 
that the WMAP data is considerably less noisy than the 
COBE-DMR dataset, and that noise can at most hide a non- 
Gaussian signal. However one should never forget the issue 
of systematics. Even though documented COBE-DMR sys- 
tematics were shown not to correlate with the observed non- 
Gaussian signal, the new data allows us to cast a new look 
at the problem. We find that a highly non-subtle combina- 
tion of QuadCube pixelization systematics and the "eclipse 
data" are to be blamed for the observed effect. 

This, we believe, closes the story. The moral is clear: 
care must be exercised regarding similar claims currently 
being made with WMAP. 



2 THE BISPECTRUM 

We start by reviewing some results and definitions pertain- 



ment high noise levels could at most dilute the observed 
signal, in spite of the noise correlations and anisotropy 
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ing to the bispectrum. Given a full-sky map, -^-(n) 
may be expanded into Spherical Harmonic functions: 
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algorithm). The most general cubic invariant is the 
bispectrum, and is given by 
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where the (. . .) is the Wigner 3J symbol. The proportional- 
ity constant is usually chosen in order to enforce a ro ughly 
consta nt cosmic variance. In Ferreira. Maeueiio fe Gorskil 
(1998) the choice was made £1 = £2 = £3, leading to the 
"single-^" bispectrum Bi = Beee- Other bispectrum compo- 
nents are sensitive to correlations between different scales. 
Selection rules require that £i+£2 +£3 be even. The simplest 
chain of correlators is therefore At — Bt- ne+i - the "inter - 
£ bispectrum - and this was studied in iMagueiid 
Other compon ents, involving more dista nt multipoles, may 
be considered dSandvikfc MagueiidEool but they are very 
likely to be dominated by noise; it is natural to assume that 
possible non-Gaussian inter-scale correlations decay with £ 
separation. 

We shall therefore consider ratios 
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where Ct — 2I+1 Em l a *m| 2 - These quantities are dimen- 
sionless, and therefore less dependent upon the power spec- 
trum. They are also invariant under rotations and parity. 

The theoretical importance of the bispectrum as a non- 
Gaussian quali f ier ha s been recognize d in a number of pub- 
lications jLuol lll994h. iPeebled (Il998l). ISpergel fc Goldberg! 
|l999 [) . lGoldberg fc Spergell Jl999D . IWang fe Kamionkowskil 
( 2000)). |Kogut_et_aj| Jl99jj) measured the pseudocollapsed 
and equilateral three point function of the DMR four year 
data. The bispectrum may be regarded as the Fourier space 
counterpart of the three point function. 



3 BISPECTRUM ANALYSIS OF WMAP DATA 

The WMAP mission jBennett fe alll2003l) was designed to 
make full sky CMB maps with unprecedented accuracy. 
There are ten differencing assemblies (DAs) in total, four in 
the W band at 94 GHz, two V band at 61 GHz, two Q band 
at 41 GHz, one Ka band at 33 GHz, and one K band at 23 
GHz. The K and Ka bands are dominated by galactic emis- 
sion and therefore neglected for cosmological analysis. The 
maps are made using the HEALPix 1 format with n side=512 
jGorski. Hivon feWandeltlll998l ; iGorski fc allll999D . The to- 
tal number of pixels in each map is 12 x nside 2 = 3, 145, 728. 
We use the coadded sum map of the Q, V and W maps, 



The HEALPix website is http://www.cso.org/sciencc/healpix 
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Figure 1. The COBE inter-£ non-Gaussian signal (bottom), us- 
ing three different data renditions: HEALPix (solid line) and 
Quad-Cube ecliptic (dashed) and galactic (dotted). In the top 
pannel we have re-reproduced the COBE galactic Quad-cube data 
(dotted), and superposed a COBE map in the same pixelization 
but without the eclipse data (dashed), and the WMAP results 
(solid). The continuum dotted lines represent 1, 2 and 3 sigma 
contours for the bispectrum expected from Gaussian maps. 
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where Ti is the sky map for the DA i with the foreground 
galactic signal subtracted, and cr^, is the no ise per obser- 
vation for DA i, whose values are given by iBennett fc all 
(2003). We use the publicly available 'foreground cleaned' 
maps, where the Galactic foreground signal, consisting of 
synchrotron, free-free, and dust emission, was removed using 
th e 3-band, 5-paramete r template fitting method described 
in lBennett et al.l i2003T) . We then use the KpO mask to cut 
the Galactic p l ane em ission and the known point sources 
iBennett et alJ l l200aft ), retaining 76.8% of the sky. 

The monopole and dipole are removed and we perform a 
harmonic analysis of the map obtaining the ai m up to / = 20. 
This is performed using the FORTRAN utility ANAFAST, 
available in the HEALPix package. We then evaluate (|1J 
and the jf obtained for the WMAP data are compared to 
the distributions P{Jf) obtained from Gaussian simulations 
subject to the appropriate beam, galactic mask and noise. 
These simulations take as input the LCDM power-law pri- 
mord ial power spectrum fit to the WMAP, CBI and ACBAR 
data dBennett et alJ ll2003f) ). The random fields were gener- 
ated using the utility SYNFAST of the HEALPix package. 
The distributions so obtained do not vary significantly from 
those obtained previously for COBE. 



4 RESULTS AND DISCUSSION 

In Fig. we plot the inter-^ bispectrum of several ren- 
ditions of the COBE data, and of WMAP. The original 
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Figure 2. The inter-£ spectrum for the various WMAP Q, V, 
and W bands (dotted lines). The solid line on all panncls is the 
spectrum for the co-added maps. 

COBE-DMR non-Gaussian inter-^ signal jMaeueiiol bood) 
was found using QuadCube pixelization. Possible deficien- 
cies of this scheme were evaluated by aligning the pixeliza- 
tion system with galactic and ecliptic coordinates. The Jf 
in both frames are plotted in the bottom panel of Fig. 
The dotted contours represent 1, 2 and 3 sigma lines for the 
bispectrum arising from Gaussian realizations. 

One would expect to see one in three coefficients ly- 
ing outside the 1-sigma contour. Instead, for COBE-DMR 
maps we see a very close alignment with the peak of the 
distribution. This gives a reduced chi squared X 2 = 0.14 
and X 2 = 0.22 for data in galactic and ecliptic pixeliza- 
tion, respectively. Computing the distributions P(X 2 ) leads 
to P(X 2 > 0.14) = 0.998 (and P(X 2 > 0.22) = 0.985) for 
maps in galactic (ecliptic) pixelization. These are the confi- 
dence levels for rejecting Gaussianity on the grounds of the 
COBE bispectrum. 

As can be seen from the upper panel in Fig. the 
WMAP bispectrum does not have such an obvious lack of 
scatter. Indeed it leads to X 2 — 0.59, consistent with Gaus- 
sianity. 

What can be the origin of this discrepancy? It is at once 
clear that we are not comparing like with like. The WMAP 
project used the HEALPix pixelization scheme, and the tests 
made for the impact of QuadCube upon the COBE map may 
not be conclusive. For example there may be an isotropic 
systematic effect, present no matter how one orients the co- 
ordinate system. To address this in Fig. Q we plotted (solid 
line in bottom panel) the jf from COBE 4 year maps ren- 
dered in HEALPix. We found X 2 = 0.26, reducing the con- 
fidence level to 96%. Hence the COBE/ WMAP discrepancy 
can partly be blamed on a poor pixelization scheme, but 
this is not enough: even in the HEALPix rendition COBE 
is much more non-Gaussian than WMAP. 

Could we have missed the non-Gaussian signal in 
WMAP by looking at the wrong combination of frequen- 
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Figure 3. The bottom pannel shows the J? for WMAP subject 
to the COBE mask (solid line) and to the KpO mask (dotted). The 
middle panel shows the J? for the COBE maps in HEALPix with 
(dashed) and without (solid) the eclipse data. The top panel (solid 
line) shows the inter-f spectrum of a map obtained by subtracting 
WMAP to COBE (both with the COBE mask). We have added 
the COBE (HEALPix) bispectrum for reference. 



cies? Equation © favours the Q band, which is the most 
contaminated by galactic emission. Galactic emissions have 
been proved to degrade non-Gaussian large angle signals 
iMagueiioll200dl . However, as can be seen in Fig. [5] there is 
a remarkable consistency between the bispectrum observed 
in all WMAP channels, a tribute to the very high signal 
to noise, but also to the lack of galactic contamination. By 
way of contrast the COBE inter-i! non-Gaussian signal came 
mainly from the 53 GHz channel, which was also the least 
noisy channel. 

A related possible source of discrepancy is the galac- 
tic mask. The KpO mask is significan tly smaller than th e 
extended galactic cut used in COBE (iBandav fc allll997l) . 
Could the extra regions near the galactic plane hide a non- 
Gaussian signal? As the bottom pannel in Fig. [3] shows this 
is not the case. Subjecting the WMAP data to the extended 
cut used by the COBE team in fact increases the WMAP 
chi squared to X 2 = 0.64. 

It would therefore appear that nothing has been missed, 
and that the WMAP bispectrum on "COBE" large angular- 
scales is indeed consistent with Gaussianity. Given that the 
higher noise levels in the COBE maps c an at most dilute 
a non-Gaussian signal (a fact proved in iMagueiiol (2000) 
even after taking noise correlations into account), we may 
conclude that a systematic error is behind the COBE non- 
Gaussian signal. 

Nevertheless, identifying the culprit is far from ob- 
vious: no documented COBE systematic mimics the ob- 
served inter-l signal. The effe cts of the eclipse data 
jBandav. Zaroubi fc Gorskl [l999 ) on the Jf, for example, 
are assessed in the top panel of Fig. They give X 2 = 0.18. 
leading a confidence level for rejecting Gaussianity of 99.2%. 
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Figure 4. The vertical thick dashed line is the value of the ob- 
served J ; 3 . The solid line histogram is the pdf of the J? obtained 
from 4020 Gaussian simulations of the sky with noise and the 
KpO mask. The dashed line histogram superimposed is the same 
pdf obtained with the COBE's extended galactic cut and DMR 



It turns out, however, that if we reject data collected dur- 
ing the eclipse season and construct a DMR 4-year map in 
HEALPix the puzzle is solved. This is shown in the middle 
panel in Fig. [31 where we have plotted the jf inferred from 
such (co-added) maps. Although the difference might look 
subtle, there is quite a substantial difference around I = 16. 
In fact the reduced chi squared in now X 2 — 0.42, consistent 
with Gaussianity. 

Curiously the single- 1 n on- Ga ussian signal found by 
iFerreira. Maeueiio fc Gorskil (I1998T) results mainly from if 
at t — 16. This stops being a severe deviant once the eclipse 
data is excluded. It would now appear that the same hap- 
pens for the Jf, but only after a better pixelization scheme 
is introduced. 



5 CONCLUSIONS 

The conclusion is that a combination of pixelization ef- 
fects and the "eclipse" systematic is behind the inter-^ non- 
Gaussian effect previously reported for COBE maps. Curi- 
ously if we take a difference map (COBE minus WMAP) the 
anomalous alignment of the jf is not present (see top panel 
of Fig. [3J . On the scales we are considering such maps pic- 
ture COBE noise, plus COBE systematics, minus WMAP 
systematics (assumed to be small). Thus one would expect 
the COBE jf signal to be enhanced in the difference map. 
The fact that it is not results from the non-linearity of the 
statistic being used, plus the subtle interplay of signal and 
systematic via the QuadCube pixelization scheme. 

In this paper we have concentrated on jf, because the 
single scale bispectrum anomalies have been explained long 
ago. However we have checked that no new anomalies on 
large angular scales emerge in the WMAP data. 



We reserve to a future publication a complete study of 
the bispectrum of WMAP of smaller a ngular scales. To our 
mind the work of Komatsu & al ( 2003) is just the beginning. 
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